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The confinement of the heavy quark 

In the previous talk [1] we have considered some aspects of the theory of 
supercharged nuclei, and came to the conclusion, that the superbound atoms 
are stable mainly due to the Pauli principle. 

Before going to the heavy quark, let us discuss briefly, how the problem 
of the supercharged nucleus can be handled practically. One possibility to 
formalize it is the following. We have to calculate the Green function of the 
Dirac equation in external field : 



V G(x,x') = -S(x-x') , (1) 

i 



where 



V = 7m (fy - iA») . (2) 

The initial external field is that of the considered nucleus. However, in the 
presence of a charge Z the stationary states correspond to atomic states with 
a charge Z — N where N is large enough to fulfil Z — N < Z cr , i.e. there have 
to be N electrons rotating around the nucleus. But if so, we will have to take 
into account, that the field which acts on each electron is also changing. 

The ground state in our case is an atomic type state. This means, that 
although we don't know what our field is, we can find it. The gradient 
squared of the zeroth component A = A of the potential has to be equal 



V 2 A = el (pnud + Pel) 



(3) 
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where p e £ =< ^70^ >• 

The charge density will depend on the potential, and the equation be- 
comes a Thomas-Fermi type equation for the potential existing in this sys- 
tem. If we find a wave function satisfying the equation we will know A, the 
Green function and p. This means, that the formal way of solution is just to 
solve a self-consistent, Thomas-Fermi type equation for the effective atomic 
potential. 

Let us consider the total charge Q 




the density is the average value 

p =< tfj( r )ip(r) >= ^ ^n(r)^ n (r) (5) 

which in the Dirac picture is the sum of all negative energy levels. The 
total charge Q will be the sum of the energy levels over all uj n . The sum 
is divergent, and we have to make a cut-off and subtract the bare particles. 
This, however, will not be enough. Indeed, if we just subtract the value of the 
charge of the free vacuum, then at Z ^ the charge of the vacuum (i.e. of 
the nucleus) will not be equal Z and will continue to diverge logarithmically. 
We have to subtract the value at small Z in a way which gives zero for 
the total charge of the electron vacuum (this corresponds to the correct 
renormalization of the charge of the nucleus). 

The above procedure does not reflect literally the subtraction of the 
vacuum charge without external field. This becomes especially obvious, if 
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one makes use of the Levinson theorem which connects the number of states 
of a particle in external field in a given energy interval with the phase of 
scattering of this particle in the same field. The number N of additional 
states in the energy interval between E\ and E 2 is defined by the difference 
of phase shift 

N = — [5{E 2 ) - 5(E 1 )] . (6) 

7T 

Because of this, the number of new states in the Dirac vacuum equals 

N = - [6(-m) - S(-oo)] . (7) 

TV 

For Z < Z c this number should be zero. However, for the Dirac equation in 
external field this condition is not fulfilled, because S(— oo) 7^ 0, and in order 
to obtain the proper definition of the vacuum charge we have to subtract a 
quantity which, generally speaking, depends on Z (S(— 00)). This subtraction 
means, in fact, that we have to change the interaction with the external field 
when E — > —00 so that S(— 00) = 0. (S(—m) can always be considered to be 
zero if the field is small). If, however, Z will be increased and becomes more 
than critical, then, as we saw, the levels will pass the point E = — m and 
move to the complex plane. It is easy to check, that every time this happens 
the phase S(—m) is changing by it and the move of n levels into the complex 
plane changes the number of states by n so that the charge of the vacuum 
(i.e. of the atom) becomes Z — n. The value of p e which enters the equation 
for the self- consistent field is to be defined by the contribution of levels which 
passed through E = —m. 

We discussed in detail this concrete structure in order to refer to it in the 
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following, talking about the heavy quark in the vacuum of the light quark. 
We shall suppose, that due to gluonic vacuum polarization the effective cou- 
pling a(r) which at small r has the usual perturbative behaviour, reaches a 
constant value at r >> ro (ro ~ 1/A), and this constant value will be more 
than unity. (Without this ansatz, allowing that a continues to grow, things 
are more complicated, but nothing essentially changes). 



oc(r) 




> r 



Figure 1 

The supposed behaviour of the effective coupling. 

It is natural to expect such a behaviour of the charge in QCD. But even an 
Abelian theory can reveal such a behaviour of charge, if it originates from a 
non-Abelian theory via spontaneous colour symmetry breaking. In this case 
the charge will increase, and as a result, 6-7 gluons acquire masses. After that 
there remain one or two massless objects - "photons" - and the behaviour of 
charge will be exactly as discussed, because after the symmetry breaking we 
have ro = ^ (heavy gluon). We can ask now : what happens in the vacuum 



6 

of light quarks under these circumstances ? Outside the region where a is 
growing, we will have an Abelian theory and we can consider the quark states 
in the normal way. 



a(r) 




Vm 
(symmetry 
breaking) 

If a becomes more than critical, the corresponding level goes to the Dirac 
sea, which, consequently, will have to be filled up and we will find an atomic 
type state. If this is so, the charge density will be exactly the same as in the 
case of a supercharged nucleus. The Xqcd is analogous to the radius of the 
nucleus. There is, however, an important difference between QED and this 
case. Indeed, in QED we consider a nucleus with a charge Z in the centre, 
and we put one or two electrons around it to organize an atom. In QCD we 
have a heavy quark with a very small intrinsic charge. 
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:_qi- 




Due to vacuum polarization the charge becomes large inside 1/A. The system 
creates an empty level, and we have to fill it. This means, that we have to add 
an intrinsic charge, equal to that of the heavy quark, but with the opposite 
sign (outside the region of 1/A). So we have here two intrinsic charges with 
the total charge zero, and of course the vacuum polarization can never change 
the total charge. In other words, in QCD the supercritical atom would be a 
meson-type state with zero colour charge. Our task is to show, that these are 
not only words - we have to include the mechanism of vacuum polarization 
formally. In order to do so, let us discuss the problem in a language similar to 
QED, neglecting non-Abelian fluctuations of the colour field A - the average 
field of heavy and light quarks inside the meson. In this case A is defined 
by the equation 

V 2 A = e 2 ( Ph - p e ) , (8) 

where ph is the density of the heavy quark, and pi that of the light quark. 
Again, we can write the equation for the Green function 



(V-to) G(x,x') = -Mx) 



(9) 
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But now the problem is, how to calculate this. What means, e. g., pheavy ? 
We can not use the expression p^ = eQ<5(r), since it does not include the 
bosonic vacuum polarization, which has to be taken into account. It is well 
known, that in principle vacuum polarization means summing a diagram like 







with gluonic loops inside. The problem is, however, how to write this in 
normal space-time language. There is a good way to sum this diagram by 
writing an equation for currents. Knowing the external current and wishing 
to calculate the total current, we have to solve the equation 

JM) = + ^J 6 '(( X - x ') 2 ) U*')d*x' , (10) 

which corresponds to the summation of the diagram. We will not derive this 
equation, because it is almost obvious for the static charge which we are 
interested in. 

Let us consider a static charge, not depending on time. In this case 

P^ = Pr^ + ^ 2 J y^ P (r) , |r-r'|>e , (11) 

where 1/e is the ultraviolet cut-off, b = ^-n c — |n/. 

The solution of this equation is very simple and leads to the usual expres- 
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sion for charge renormalization in QCD. In order to see this, let us introduce 
the quantity Q(r) : 



Q(r) = [ p(r')4Tvr' 2 dr' . (12) 
Jo 

The logarithmic derivative of Q(r) is 



dzQ(r) = r^p- = 4nr 3 p(r) . (13) 
For Q(r) we can write the equation 



(14) 



The integration in the right hand-side of this expression contains two loga- 
rithmic regions : \r — r'\ << r and r >> r' . In the first region, d^'Q(r') 
can be substituted by d^Q(r) ; integrating over the second region, r' in the 
denominator \r — r'\ 3 can be neglected. As a result, we obtain 

d,Q(r) = d, Q ext {r) + g jf Q(r) + ^Q(r) , 

which is equivalent to 

d c (l - Q(r) = Q«rt(r) (15) 

or 
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Q{r) = Qext [ r) . (16) 

2n e 



For a point-like charge 



Qext(r) = 1 



and we have 



A (r) = ^ = a(r)A ext . (17) 



The concrete expression 

a {r) = ^ (18) 

Ait e 

obtained from perturbation theory has, of course, an unphysical singularity. 
For a point-like charge (18) has to be substituted by an expression corres- 
ponding to the behaviour of a{r) as shown in Fig. 1. However, for the 
distributed charge the relation between the external field and the field which 
takes into account the polarization is non-local in coordinate space. The cor- 
rect expression for the relation between the external field and the observable 
field is local in the momentum space : 

A(q) = a(q 2 )A ext (q) (19) 
which leads in the coordinate space to an expression of the following type : 

A(r) = J K(r-r')A ext (r')d 3 r' , (20) 
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where 



K(r) = j e iqr a(q) 



d 3 q 
Similarly, 



p(r) = J K(r-r')p ext (r')d 3 r' . (21) 

If we now suppose, that a(q) as a function of 1/q behaves according to Fig. 1, 
and at large q values it is defined by perturbation theory, then the equation 
(8) has to be understood as an equation for A ext : 



V 2 A ext (r) = 5(r) - ^(r) 7o ^(r) (22) 
where ipe(r) is the solution of the Dirac equation 



(H + A)tp = Eip 

for the light antiquark in a superbound state in the field A{r) defined by 
equation (20). The solution of this problem gives the energy and the features 
of the meson q n q£ with zero total charge. Due to (21), 



Q = J K(r)d 3 r Q ext = a(q = 0)Q ext (23) 

with Q = if Q ext = 0. 

We have just proved, that because of the big charge which appears 
through vacuum polarization, in the case of QCD the "atomic" bound state 
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will, indeed, be a meson. The heavy quark will decay into a g^-meson and 
a light quark : 

q h ^M (q h q £ ) + q k . (24) 

In the next lecture we shall consider light quark states. So far there is one 
important thing to stress, namely if we don't include essential interactions 
between light quarks in the vacuum, we come to a reasonable conclusion for 
the case of the heavy quark but, as we will just see, to an unreasonable one 
about the light quark. Indeed, let us try to extend the considered procedure 
to the latter case. Suppose, that there is a light quark moving, and a potential 
acts on its vacuum. What will we see classically ? Since the Coulomb field 
is a vector field, it is shrinking, but the total integral remains the same. 
Because of this, we will find immediately, that there is a bound state in this 
potential even for fastly moving particles. This, however, means, that we 
have here an unstable state, which has to be filled, and as a consequence the 
light quark will decay into a meson and a light quark again : 

qe — > M + qn 

which, of course, contradicts the energy conservation, unless the appearing 
meson is of negative energy. In order to have a self-consistent picture, we 
have to suppose that the light quark in the vacuum will interact so strongly 
that there have to be negative energy levels and the whole vacuum has to 
be rearranged. So from the picture we described we come quite naturally to 
light quark interactions. We will see, that these interactions are, indeed, very 
strong and lead to the confinement of light quarks which will take place at 
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relatively small ^ values ; this means, that the overall corrections for vacuum 
polarization will not be large. 

This is for the future. What we have to add now, is, that even in the 
language which was accepted so far, with no strong interactions between 
particles in the vacuum, the problem in real QCD which is non-Abelian is 
more complicated. In QED we have one charge, and all the electrons in the 
vacuum interact with this charge, independently from each other. In QCD 
this can take place only, if the field of the heavy quark is an Abelian one. It 
is highly probable, that this is, indeed, the case, when the field of the heavy 
quark becomes large as a result of gluonic vacuum polarization. 
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